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ABSTRACT 



We analyze the physical parameters of interstellar filaments that we describe by an idealized model of isothermal self- 
gravitating infinite cylinder in pressure equilibrium with the ambient medium. Their gravitational state is characterized 
by the ratio /cyi of their mass line density to the maximum possible value for a cylinder in a vacuum. Equilibrium 
solutions exist only for f^yi < 1. This ratio is used in providing analytical expressions for the central density, the 
radius, the profile of the column density, the column density through the cloud centre, and the fwhm. The dependence 
of the physical properties on external pressure and temperature is discussed and directly compared to the case of 
pressure-confined isothermal self-gravitating spheres. Comparison with recent observations of the fwhm and the central 
column density Ah(0) show good agreement and suggest a filament temperature of ~ 10 K and an external pressure 
Pext/k in the range 1.5 x 10* K cm"'' to 5 x 10* K cm~^. Stability considerations indicate that interstellar filaments 
become increasingly gravitationally unstable with mass line ratio /cyi approaching unity. For intermediate /cyi > 0.5 
the instabilities should promote core formation through compression, with a separation of about five times the fwhm. 
We discuss the nature of filaments with high mass line densities and their relevance to gravitational fragmentation and 
star formation. 
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1. Introduction 

Filamentary structures are an ubiquitous phenomenon 
in the interstellar medium. Thanks to the high an- 
gular resolution and signal to noise of dust imaging 
with the Herschel Space Observatory, it is now possi- 
ble to quantify the basic empirical properties of fila- 
ments. Filaments have been imaged in the submillime- 
tre in exquisite detail in non star- forming such as Polari s 
(jMen'shchikov et al.l I20T0I: iMiville-DeschSnes et al1l20l0ll . 
in molecular regions with low-mass st ar formation such as 
Aqui l a and IC 5146 ( And r e et a l..,2010i;lArzoumanian et al.l 



2010 ) are observed mainly along filaments (lAndre et a.1.1 



20101: IMen'shchikov et all 120101: lArzoumanian et al.l 1201 iD 



I2OIII: IMen'shchikov et al.l 12010( 1. and in higher-mass re 
gions including Vela C and the Rosette molecular cloud 
()Hill et al.ir2011l:ISchneider et al.ll2012[) . At least in the low- 
mass star forming regions and i n Polaris , the characterist ic 
filament sizes are about 0.1 pc (jArzoumanian "eraI1[20TTh Fl 
Clear evidence for star formation seems to be re- 
lated to cloud structures with a typical extinction larger 
then Ay ^ 7 ma g (Onishi ct al. 1998; Johnstone et al.l 
12004 lAndre et al.l feOlO'). which corresponds to a col- 
umn density A'h about 0.6 x 10^^ cm~^ (see com- 



ments in Appendix IA)) n Furthermore, the Herschel stud- 
ies show that high column densities are typically as- 
sociated with filaments and, finally, that cold dense 
clouds called "prestellar cores" possibly related to the 
very early stages of the star formation (jKonvves et al.l 



^ Appendix |X| discusses uncertainties in this and other impor- 
tant physical parameters like column density. 

^ According to lEnoch et al.l (|2007f) . however, such a threshold 
might vary from region to region. 



Empirically, it is found that filaments that are star forming 
are characterized by not only a higher central colu mn den- 
sity but also a higher ma s s per unit length ()Andre et al.l 
I2OIOI : lArzoumanian et aD 120111 ). It is argued that the 
cores are possibly created out of the elongated structures 
through gravitational in s tabilities ^S chneider & Elmcgrce 
1979; Gaida ct al. '1984; Han awa et al I il993t ICurrvl 120^ 
A ndre ct al. 2010; ,AIen'shchikov et al.ll2oTo[ l 

It is important to develop a physical model for the 
filaments so that the new measurements can be exploited. 
That is the goal of this paper. The basic model that 
we explore is an isothermal cylinder confined at a finite 
boundary rcyi by an external pressure Poxt provided by 
an ambient or embedding medium. This is certainly 
more realistic than an isothermal cylinder in vacuum, 
whose properties (such as the radial profile) are not 
surprisingly in disagreement with observations (e.g., 
IMen'shchikov et al.l [20101: lArzoumanian et~aI1l201ll ). Even 
these pressure-confined models are unlikely to describe 
what is a very complex interstellar medium where the 
origin of the filamentary structures is still a mystery, but 
the models can serve the same role in developing our 
understanding that has been played by the spher ically- 
symmctric analog, the Bonnor-E b ert sphere (lEbert 
1955; Bonnor 1956; Nag asawa [1987 : Inutsuka fc Mivama 
1992, 1997; Curry & McKee* '2000'; 'Fiege & Pudritj 
2000,; Kandori_et al .2005; Fischera fc D opita 200.^. 
Furthermore, we find that the models do actually provide 
a consistent description of the sizes and column densities 
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recently reported for filaments and some insight into the 
instabilities relating to filaments. 

We analyze systematically the physical properties of 
pressure-confined filaments, comparing and contrasting 
these results with the corresponding properties for spheres. 
We develop analytical solutions for several observables. 
Parts of the results, on which we build, can be found in 
earlier work o n the spectral energy di stribution of inter- 
stellar clouds (jFischera k, Dopital l2008l - paper I) and on 
the spectral energy distribution of condensed core s embed- 
ded i n pressurized elongated or spherical clouds ()Fischeral 
donl- paper II). 

Our paper is organized as follows. The physical model 
is presented in Section [51 from which properties are derived 
in Section [21 Observed properties corresponding to size and 
column density are compared with the model in Section [31 
with remarkable agreement. Structure formation due to 
instability of the filaments is discussed in Section [5l It 
seems plausible that dominant protostellar cores arranged 
alon g filaments could arise from a "compressive instability" 
(jNag asawa 1987) in high column density filaments. Some 
other instabilities seem less relevant. A summary and dis- 
cussion is provided in Section [SI Several appendices deal 
with the astronomical context and details of some calcula- 
tions. 



2. Physical model 

We assume that the filaments can be described by a model 
of isothermal self-gravitating infinitely long cylinders. For 
the isothermal gas, the relationship between gas pressure p 
and gas density p is 



K = p/p = Tk/{pmii) 



(1) 



where k is the Boltzmann constant, /i the mean molecu- 
lar weight (2.36 for a molecular cloud with cosmic abun- 
dances), and toh the mass of the hydrogen atom; the 
isothermal temperature T is considered to be an effective 
value related to both the thermal and the turbulent motion 
of the gas (see Appendix IA.4[) . If the effective temperature 
is equal to the thermal or kinetic temperature then K = c^, 
the square of the sound speed. Furthermore, as discussed in 
the introduction, the clouds are assumed to be in pressure 
equilibrium with the surrounding medium of pressure Pcxt ■ 
This does not imply that the external medium is isother- 
mal or at the same temperature; it is simply a boundary 
condition. 

2.1. Cloud profile 

Isothermal self-gravitating infinitely long cylinders have a 
well-known density o r pressure profi le given by (see e.g. 
IStodolkiewic^ (ll963h : IOstriked 



p{r) 



where 



(1 + (r/V8ro)2))2 



(2) 



rl = Kl{A^Gp,), 



(3) 

Pc is the central density, and G the gravitational constant. 
For pressure-confined clouds the profile terminates where 
the cloud pressure is equal to the external pressure Pcxt] this 
defines the cylinder radius Tcyi. As in papers I and II, in the 
following we refer to the ratio Pc/pcxt as the "overpressure." 



2.2. Mass line density 

Analogous to the cloud mass in the case of spherical clouds 
we have in the case of cylinders the mass per cloud length 
(mass line density), given by 



M 



dr 2'Kr p{r) 





2K 



l + reyiV(8r2) 



(4) 



In the limit Tcyi 3> V^tq the mass per cloud lengths ap- 
proaches asymptotically a maximum value given by 



2K 
~G 



= 16.4 



T 
10 K 



Mo pe- 



ls) 



For pressurized clouds this gravitational state corresponds 
to an infinite overpressure, but there is no dependence on 

Pcxt- 

For comparison, in the case of an isothermal self- 
gravitating pressurized sphere, the maximum mass is iden- 
tical to its critical mass, given by (see Appendix [Bjl 



sph.cr 



4.191- 



V47rG3poxt 



(6) 



This gravitational state is characterized by an overpressure 
Pc/Pcxt ~ 14.04. As shown in paper II, a critical stable state 
is related to the maximum possible pressure that a sphere 
of given mass and K (temperature) can produce at the 
cloud outskirts under compression. If the external pressure 
is increased slightly, the resulting equilibrium configuration 
is unstable to gravitational collapse, and so this is aptly 
called a "critical stable sphere" . The spherical clouds can be 
characterized as "subcritical" or "supercritical" depending 
on whether the overpressure is lower or higher than the 
value for a critical stable sphere. Note that for overpressures 
higher than 14.04, the mass of an equilibrium configuration 
is smaller than the critical mass. Thus we can also say that 
for given pcxt and K the critical mass is the maximum mass 
for which an equilibrium solution exists. 

Unlike Msph,cr, the maximum mass line density depends 
only linearly (compared to quadratically) on K or cloud 
temperature and does not depend at all on the size of the 
external pressure confining the cloud. We discuss the basic 
consequences in Sect. 13.61 and Sect. 13.71 

There are no equilibrium solutions for cylinders or 
sphe res with masse s above the corresponding maximum val- 
ues (|McCrealll957[ ). They would collapse radially towards 
a line (spindle) or to a singular point, respectively. This 
has led to the terminology "supercritical mass" when M > 
-Msph,cr, or "supercritical filament" when AI/l > (M//)max- 
The gravitational states considered in this paper are all 
equilibrium solutions. We prefer to reserve the term "su- 
percritical" to describe the equilibrium states for spheres, 
as above, and note the fundamental difference that there 
are no such supercritical equilibrium states for cylinders 
(Sect.[ 



2.2.1. Mass fraction fcy\ 

A very useful quantity is the mass fraction 



/cyl - {M/l)/{M/l\ 



< 1. 



(7) 
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For cylinders wc find directly from Eqs. [T] to [7] that the 
overpressure 

PcMxt=l/(l-/cyl)'. (8) 

We can therefore use /cyi to characterize the gravitational 
state of the cylinders. This mass ratio therefore appears in 
the formulae below for the derived properties of cylinders in 
different gravitational states. Two instructive limits to con- 
sider are the case of high overpressure or /cyi 1, which 
we shall see corresponds to a narrow high column density 
filament with strong self-gravity, and the case of vanishing 
overpessure or /cyi — > 0, corresponding to a filament much 
more difficult to recognize or characterize. There is, how- 
ever, no distinction between "pressure-confined" and "self- 
gravitating" cylinders, since both features are intrinsic to 
all equilibrium models. 

This parallels the approach used in paper I where the 
gravitational state of pressurized spheres was described by 
the mass ratio /sph = M/Msph,ci- of the cloud mass relative 
to the mass of a critical stable sphere. While in paper I 
the physical parameters for spheres were discussed only for 
the physically stable regime with overpressures below the 
critical value 14.04, we consider here the whole physical 
range. This leads to the curious situation that for large 
overpressures above the critical value the gravitational state 
is described by /sph < 1. 

3. Derived properties 

In the following wc describe the physical properties of 
molecular cylinders and compare spheres for a given cloud 
temperature and external pressure. We adopt illustra- 
tive reference values for temperature and pressure. See 
Appendix |2] for further discussion. For the cold molec- 
ular clouds we adopt K corresponding to 10 K (cs = 
0.19 km s~^). For the external pressure we adopt Pcxt/k — 
2 X 10^ K cm^'', consistent with paper I and paper II. 

3.1. Central density 

For the central density we have from Eqs. [T] and [8] 



10- 



Poxt 



The central number density nji{r = 0) = nH(0) is obtained 
by evaluating 

nii = P (y^^ ^fiiiriB^ =p{p.mB_y^, (10) 

where ni/nn is the relative abundance of element i com- 
pared to H, and /i^ is its relative mass units of ttt-h; on the 
right, the effective mass per H, p., is 1.4 for cosmic abun- 
dances. Figure [1] shows these central values as a function of 
the mass fraction. 

For a given temperature and external pressure the cen- 
tral density increases as 1/(1 — /cyi)^ and so in the limit 
/cyi — > 1 the density becomes infinite. 

For a given mass fraction, the density in the cloud cen- 
tre is proportional to PcxtfJ-/T. Cylinders have a higher 
overpressure or overdensity in the cloud centre compared 
to spheres (Fig. [1} for the same mass fraction, except in 
the case of small mass fractions (/ 0.15). In cylinders 
the overpressure of 14.04 corresponding to critical stable 
spheres is achieved for a mass fraction of /cyi ~ 0.733. 



10' 



10' 



: Pex,/*:=2x10''K/cm' \ 


■ 

cylinder 




^ - 
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Fig. 1. Central density of an isothermal self-gravitating 
molecular cylinder (black) and sphere (grey), for fixed val- 
ues of the ambient pressure. Pcxt/k — 2 x 10** K cm""^ 
and cloud temperature T = 10 K. The dotted line is the 
value achieved by a critical stable sphere. This same cen- 
tral density (or overpressure) is reached in a cylinder with 
/cyi — 0.733 (filled black circle). 

3.2. Cloud radius 

Solving Eq. |4]for r^, using Eq. [9]and recasting Eq. |3]as 

(l-/cyl), (11) 



K 



ro 



provides 



r-cyi 



\/47rGpc 



/cyi 2 



(12) 



This along with the boundary radius of the pressure- 
confined sphere is plotted as rd in Figure [5] as a function 
of mass fraction. 




Fig. 2. Like Fig.[Tl but for the radius of an isothermal self- 
gravitating cylinder (black) and sphere (grey). 



The size of the cloud for given mass fraction / is propor- 
tional to (T/ fi) I Y^poxt ■ The radius for cylinders is symmet- 
rical around / = 0.5 where the cloud shows its maximum 
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size. The filament grows in size through accretion of new 
material as long as the mass ratio /cyi < 0.5. Further ac- 
cretion will cause a shrinking of the filament. In the limit 
/cyi — > 1 (the mass-line density approaches the maximum 
value 2K/G) the cylinder will become infinitely thin. 

Spheres by contrast have a finite size at large mass ra- 
tios. As pointed out in paper I the spherical cloud has a 
maximum size at /sph ~ 0.870. Pushed above this value 
by accretion an isothermal cloud size shrinks to the size 
of a critical stable cloud. Supercritical stable spheres have 
smaller sizes still. 



3.3. Radial profile of the pressure or density 




- r I r,.. 



Fig. 3. Pressure profiles of isothermal self-gravitating 
cylinders (black) and spheres (grey), for various overpres- 
sures Pc/Pcxt from 2 to 20. The dashed curves show den- 
sity profiles for simple power laws normalized at the cloud 
boundary. 



The radial dependence of pressure given by Eq. [5] is 
shown in Fig. [3] for various overpressures, which according 
to Eq. [HI are equivalent to 1/(1 — /cyi)^- For comparison, 
profiles are also shown for spheres for the same set of over- 
pressures; in the central region of the cloud their density 
falls off more strongly towards larger radii. 

In the limit of r"^ /{Sr^) ^ 2 the clouds show a steep den- 
sity profile at the cloud edges with p{r) oc r~^. However, 
this property requires clouds with very high overpressures 
(large /cyi, close to the maximum possible mass-line den- 
sity) . n the outer half of the cylinders the density profiles for 
overpressures less than 20 are shallower than a profile. 
For overpressures in the regime 6 to 12 the density pro- 
file at the outer half is more consistent with a r~^-profile. 
Compared to spheres, 

A useful characteristic size for describing properties of 
the cylinder is Vcyi for small /cyi and \/8ro for large /cyi. 



3.4. Projection of the density profile 
3.4.1. Column density profile 

The column density at an impact parameter x (in units of 
cloud radius rd) is given by the integral 



7Vh(x) = 2 



nii{0) dr 



(l + a;2r2yi/8r2 + r2/8r2)2 
The integral evaluates to: 

c L + c 



1 u 
— arctan — = 
/c Jc 



}■ 



(13) 



(14) 



with c = l+a;^''cvi/8ro and u = -v/T^-^J^cyi/ (\/8?'o)- Using 
the replacement r-^yi/Sr^ — /cyi/(l — /cyi) it is straightfor- 
ward to show that the column density profile for given /cyi 
is given by: 



Pcxt 



V8 



inG 1 - /cyi 

1 - /cyi 
1 - /cyi + a;2/cyi 



v//cyl(l-/cyl)(l-a:2) 




(15) 



The profile shape is solely determined by the mass frac- 
tion /cyi. The amplitude of the profile with fixed /cyi is 
proportional to ^ypcxt and is independent of the cloud tem- 
perature. 

In the limit /cyi — > 1 the profile for x <C 1 can be ap- 
proximated by: 



Pcxt 



V8 



47tG 1 - /cyi 

/ 1 - /cyi 

{l^f,y,il-x^) 
In the limit of low /cyi 



3/2 



(16) 



Nnix) « (pext/^)(/ir7iH)-^2Vr^rcyi, (17) 
as for a uniform-density cylinder. 

3.4.2. Central column density 

An important characteristic of self-gravitating clouds is the 
column density through the cloud centre, A^h(O), which for 
cylinders is a trivial solution of Eq. [Tn]for a; = 0: 



A^h(O) 



Pcxt 



V8 



AttG 1 - /cyi 



(/imn) 



y^/cyi(l - /cyi) + arctan ^ 



/cyi 



-/ 



cyi 



(18) 



Fig. m shows that the column density increases more 
strongly with mass fraction in the case of cylinders. 
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Fig. 4. Like Fig. [Tl but for the central column density (left panel) and mean column density (right panel) through an 
isothermal self-gravitating cylinder (black) and sphere (grey), as a function of mass fraction /. The dashed line in the 
left-hand plot is the approximation given in Eq. 1191 The dashed lines in the right-hand plot echo the central column 
densities from the left-hand plot for comparison. The horizontal dotted lines give the values achieved by a critical stable 
sphere. The values for a cylinder with that overpressure, mass fraction 0.733, are shown as filled black circles. 



In the limit of high mass fraction {fcyi >~ 0.6) the 
column density can be approximated by 



^h(O) 



Pcxt 



47rG 1 - /cyi 



^(m)-^ J = nH(0)y8ro J. (19) 



For given external pressure the central column density in- 
creases as 1/(1 — /cyi)- In the limit of maximum mass frac- 
tion the column density would be infinite. If we replace the 
mass fraction through the overpressure we find that the 
column density for high overpressure increases, as found in 
paper II, proportional to The limit for low /cyi follows 
from Eq. [T71 

Column densities in the molecular environments in 
question are measured via near-infrared colour excess or, 
in the case of the filaments being discussed, the submil- 
limetre optical depth (Appendix Nevertheless, Ay is 
widely used as "shorthand" for column density, which we 
do here as well, even though it is not a direct observ- 
able in dense molecular clouds. A column density of H 
nucleons of 10^^ cm~^ is taken to correspond to Ay — 
5.3 mag f Appendix [A)) FI In other units, this is equivalent 
to 2.3x 10~^ gm cm~^ or l.lx lO^M© pc~^. This conversion 
relates the left and right vertical scales in Fig. HI 

For the adopted external pressure, the column density 
lies between 2 x 10^^ and 2 x 10^^ cm~^ {Ay in the range 1 
to 10 mag) for mass fractions /cyi in the range roughly 0.2 
to 0.8. 

Cylinders with the same overpressure pc/Pext as spheres 
are characterized by a lower central extinction (paper II). 
For example, for the adopted external pressure, the central 
extinction through a cylinder with the same overpressure 
as a critical stable sphere is approximately Ay ~ 6.4 mag 
compared to 8 mag for that sphere (Fig. U]). 



3.4.3. Average column density 

For unresolved filaments we also give the mean column den- 
sity (see also paper II) . The mean value for cylinders is given 
by 



(tVh) = im) 



-1 M/l 



2r, 



cyi 



(jUTOh) ^ 



47rpcx 
8G 



/cyi 



1-/. 



cyi 



(20) 



We see that the mean column density of cylinders also ap- 
proaches infinity asymptotically in the limit of large mass 
ratio. 

By contrast, the mean column density even for a super- 
critical sphere is finite. In the limit of large mass fraction 
the mean value is given given by 

M 



' sph 



47rpcx 

TT^G 



■V2. 



(21) 



-'Vh2/(10^^ cm ^) is therefore approximately A\ 



As shown in the App. [Bl the value is exact where a cloud 
of given mass produces a pressure maximum at the clouds 
outskirts. The value is for example exact for the special 
case of a critical stable sphere or in the limit of infi- 
nite overpressure. For Pcxt/k = 2 x 10^ K cm~^ we have 
(iVn) « 4.36 X 10^1 cm'^ {Ay « 2.33 mag). As shown in 
Fig. 13] stable spherical clouds above /sph > 0.1 have mean 
column densities not smaller than 10^^ cm~^ and all spher- 
ical clouds, even if we include the supercritical cases, mean 
column densities not larger than 6 x 10^^ cm~^. 

3.4.4. Column density profiles for different central Ay 

According to the previous sections, for a given external 
pressure the column density or extinction in these inter- 
stellar clouds can be specified by their overpressure or /, 
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without explicit specification of cither the cloud tempera- 
ture or their mass line density (filaments) or mass (spheres). 
For filaments the extinction is related monotonically to the 
mass ratio /cyi. The same applies to spherical clouds up 
to overpressures Pc/Pcxt ~ 14.04. The extinction varies as 




Fig. 5. Column density profiles of isothermal self- 
gravitating cylinders (black lines) and spheres (grey lines) 
for several central extinction values from Ay = 1 mag 
to Ay = 8 mag. The clouds are pressurized by Pext/k = 
2 X 10"* K cm~3. 



20 



15 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
- /)exi/* = 2x10''K/cm^ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

\ ; 


sphere ^ 


"^^^\\ cylinder 




^6" \ \ 




- \\\ : 





















1.0 O.i 



0.6 0.4 0.2 0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 6. Pressure profiles (density profiles) for the column 
density profiles shown in Fig. [5l Recall that profiles for the 
same overpressure are compared in Fig. [3l 

In Fig. [5j for the adopted Pcxt/fc = 2 x 10"* K cm~'^, 
we give the profiles of the column densities of cylinders 
and spheres for a number of different values of the central 
extinction. The column density profile steepens as a func- 
tion of the central extinction or /. For the same central 
extinction, the column density of a cylinder falls off more 
quickly with impact parameter x compared to a sphere. 
This relative behaviour is opposite that for the pressure 
profile (Fig. El). 



For the same central extinction, cylinders have a higher 
overpressure Pc/Pcxt compared to spheres, as discussed in 
Sect. 13.4.21 and shown in Fig. [SI 

3.5. The fwhm 

Observationally the fwhm has more relevance than the 
physical radius which is either not directly observable or dif- 
ficult to determine because of the fluctuating background. 
The impact parameter Xf^hm of the fwhm is given by 



Nii{xf^h,n) = ^A^h(O). 



(22) 



The fwhm of cylinders and spheres is examined in 
Appendix [Cj For the adopted external pressure, values are 
shown in Fig. [71 
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Fig. 7. Like Fig. [51 but for the fwhm of an isothermal self- 
gravitating cylinder (black line) and sphere (grey line) . The 
dashed line gives the asymptotic behaviour of the fwhm of 
an isothermal self-gravitating cylinder for high / (Eg. Y2^ . 
The vertical dashed-dotted lines mark the mass fraction of 
the maxima of fwhm (Eg. I25p . The horizontal dotted line 
shows the fwhm of a critical stable sphere. The fwhm of 
a cylinder with this overpressure 14.04 (/cyi = 0.733) is 
shown as the filled black circle {fwhm^y^ = 0.0812 pc). 

For the limit /cyi — 1 an asymptotic behaviour of the 
fwhm of the cylinder can be found. Using the approximation 
of the column density profile (Eq. [T6t leads to 



Mm^yi w 2v/22/3 - iVSro. 



(23) 



which decreases linearly as (1 — /cyi)- As seen in the Fig. [3 
this behaviour is valid for /cyi >~ 0.7. If we replace the 
mass ratio by Eq.[9lwe see that the fwhm decreases inversely 
as the square root of its overpressure Pc/Pcxt- The same 
behaviour is valid for spheres (see App. [C|) . 

For the limit /cyi the asymptotic behaviour is 



'cyi' 



(24) 



which varies as y/ /cyi . 

In between these limits, the fwhm shows a maximum 
size for both spheres and cylinders, just as seen for rd. 
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Because of the steepening shape of the column density pro- 
file towards higher / the maxima appear at lower / then 
that for the maximum rd. The maxima are given by 



fwhm^ 



T 
10 K 



2 X 104 K cm 



-3 



0.145 pc, /cyi = 0.331 
^ 0.181 pc, /,ph = 0.521 



(25) 



We also note that for spheres the decrease of the fwhm from 
the maximum to that of the critical value is also stronger 
than in case of the cloud radius. The value of a critical 
stable sphere is given by: 



fwhmfpY^ — 0.106 



T 
10 K 



Poxt/fc 



2 X 104 K cm" 



pc. (26) 



3.5.1. The fwhm- Nii{0) relation 



There is a relationship between the two observables fwhm 
and iyH(O), as shown in Fig. [5] for both spheres and cyhn- 
ders. A'h(O), the column density corresponding to the above 
fwhm^^^, is given by 



7Vh(0) = 10^ 



2 X 104 K cm" 
3.524 cm-2, f^yi = 0.331 
4.552 cm-2, /^ph = 0.521 



(27) 



which is independent of K, and using Eq. (^Dfor cylinders 



fwhrn^^ w 0.051 



T 



10 K 



1 X 1022cm-2 
^h(O) 



pc. 



(28) 



For high mass ratios {fcy\ >^ 0.7) if we eliminate the 
dependence on (1 — /cyi) in approximation Eq. [231 for the 
fwhm using the inverse dependence in approximation Eg. 1191 
for the central column density we obtain the simple rela- 
tionship (see also Eq. IC.9I in Appendix [C|) 



fwhm^^i « 2v/22/3-l(i^/G)/(A^H(0) Amn), 



(29) 



which is independent of both pext and /cyi but does re- 
spond to K. For molecular gas, we have in the case of high 
overpressure 



0.111 



T 
10 K 



1 X 1022cm-2 



pc. 



(30) 



Note that this is identical to Equation[28l except for a larger 
prefactor. 

At small overpressure (flat density profile) the fwhm and 
A^h(O) increase for given temperature and external pressure 
proportionally with the cloud size. We obtain the asymp- 
totic behaviour for low mass ratio or low overpressure (see 
also Equations [HI [H and 



fwhm^yi{fcyi 0) 



V3 K 

2 Pcxt 



iVH(0)/2TOH. 



(31) 



For intermediate values of the mass fraction, the fwhm- 
A^h(O) relation for both of cylinders and spheres can be 



0.25 



0.20 



p,„/fc=2x10''K/cm' 



r,| = 10 K 
[I =2.36 



sphere 




100.0 



^h(O) [10^ 



Fig. 8. Relation of the fwhm and the central column den- 
sity A'h(O) through an isothermal self-gravitating cylinder 
(black) and sphere (grey). The cloud temperature is 10 K 
and the external pressure is 2 x 10* K cm""^, except for the 
long-dashed curve where half the pressure is assumed. The 
column densities corresponding to the maxima in the fwhm 
are marked with vertical dashed-dotted lines. The asymp- 
totic behaviour for low and high overpressure is shown by 
short-dashed curves. The dotted curves (almost hidden by 
the solid curves) are approximations using Eq. 1321 The val- 
ues for a cloud with the overpressure of a critical stable 
sphere are shown as filled circles. 



interpolated to high accuracy between the above limiting 
behaviours, as discussed in Appendix |D] The relation is 



fwhm{t) = fwhm^^^C 



(t/toT- 



(32) 



where t — Nb_{0)/ Nh{0) is the normalized central column 
density, a, b, to, and C are appropriate constants described 
in the Appendix, and the power 7 is adjusted to provide 
the best fit. The corresponding values for molecular cylin- 
ders and spheres are given in Table [1] for T = 10 K and 
the adopted external pressure (denoted with superscript 
0). For cylinders and spheres the analytical function is an 
excellent approximation of the numerical solution overplot- 
ted in Fig. |8l For cylinders the approximation also provides 
the correct power law behaviour for small and high over- 
pressure, by construction. For the spheres we allowed for 
a variation of index of the power law behaviour for large 
overpressure; however, the change is small. 



[pc] 



flO^^cm- 



7 



C 



cyi. 0.144 3.581 1.862 2 0.000 2.106 

sph. 0.181 4.648 1.796 2 -0.044 2.126 

Table 1. Derived parameters of the approximation Ea.l32] 
of the fwhm-N'ii{0) relation for spheres and cylinders. The 
temperature is 10 K and the external pressure Poxt/fc = 
2 X lO* K cm~^; parameters for adopted values are denoted 
by superscript 0. 
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From Appendix [D] follows the scaling relation of the 
column density at the maximum 



A^g(O) V2xl04Kcm-^ 
and the scaling relation for the maximum fwhm 

fwhm^^^^ \10Kj V2 X 104 K cm" 



a + b 



(33) 



(34) 



Combining these with the values in Table [T] we recover 
Equation mi 

For constant pcxt, ^h(O) is also fixed. The dependence 
of both asymptotes and fwhm^-^^^ on K or T is the same 
and so the locus of fwhm vs. A^h(O) is simply scaled by K. 

3.6. Dependence on external pressure 

In the previous sections we have examined how the phys- 
ical properties of isothermal self-gravitating cylinders and 
spheres, subject to a fixed external pressure and with a 
fixed K or temperature, change with the mass, parameter- 
ized through the mass ratio /. We can therefore directly un- 
derstand how the properties would be affected by a change 
of mass, say through accretion of additional material. 

In this and the following subsection we discuss how the 
physical parameters are affected by altering the external 
pressure and the temperature (through cooling), respec- 
tively. 

3.6.1. Changing external pressure 



As we have seen in Sect. 12. 21 the maximum mass line density 
of a cylinder does not depend on the external pressure. 
Therefore, in the case of constant mass line density /cyi and 
fixed temperature, a change of the external pressure has 
no effect on the gravitational state. In a quasi-equilibrium 
situation a cylinder can always compensate for a higher 
external pressure through compression; the central pressure 
rises so that the overpressure (Eq. |5]) remains the same. 
While a change of the external pressure has no effect on 
the density profile it does affect the appearance because a 
cylinder becomes more opaque in higher pressure regions in 
accordance with A'^h(O) oc ^pcxt fEg. [T51). 

The case of a sphere with constant temperature and 
constant mass M is quite different; a change of the external 
pressure affects the gravitational state. The ratio of M to 

the critical mass changes as /sph = \/pcxt/Pcxt : where Pc"^ 
is the maximum possible pressure that the cloud of fixed M 
can sustain at its edge (see Eq. [H]). As the external pressure 
is increased the central pressure reaches a finite maximum 
14.04pcrit when /sph = 1, after which there is no stable 
configuration. 

Furthering the discussion in Sect. 12.21 a cloud is consid- 
ered to be stable if the compression leads to a pressure in- 
crease at the cloud edge, i.e., the response —(dp(rci)/drci) > 
where is the cloud radius. Subcritical spheres, those 
with overpressure below 14.04, are stable against compres- 
sion with a pressure response — (dp(rsph)/drshp) > 0. Just 
above the critical state, a sphere is unstable|3 However, we 



also note that supercritical spheres significantly above the 
critical overpressure are not all unstable according to the 
above definition. 

In the case of a cylinder the response is simply given by 



dj'(''cyl) ^ r, P{rcy\) ^ Q 



dr. 



cyi 



rcy\ 



(35) 



for all mass ratios < /cyi < 1. Because of this, all gravi- 
tational states are stable. In the context of the equilibrium 
solutions being discussed, there is no critical state (or su- 
percritical states) for cylinders as there is for spheres (see 
also ,Ebert.l955; .McCrea.1957.1 . 

3.6.2. Embedded filaments 

Filaments are often characterized by comparing their mass 
line density to the maximum value 2K/G. A basic problem 
of interpretation which arises concerns the effective temper- 
ature of the filament, in particular if the filament itself is 
not entirely isothermal. For example, larger filaments with a 
relatively high effective temperature of 50 K or more might 
form/ contain small dense interior filaments which are a fac- 
tor five or more colder. The profile of the column density 
would show a narrow distribution on top of a broader one. 
The small size of the narrow feature is related to the colder 
gas temperature and to its higher external pressure, which 
in this case is the pressure in the interior in the broad fila- 
ment. Depending on the distribution of mass assigned, both 
broad and narrow filaments might be below their relevant 
maximum mass line density (the fact that they were ob- 
served would suggest this conclusion). See further discus- 
sion in Sects. 14.21 and 15.31 



3. 7. Dependence on temperature 

To understand how the isothermal self-gravitating clouds 
are affected though cloud cooling we consider a spherical 
cloud of fixed mass and a cylindrical filament with fixed 
mass line density. For a given external pressure, equilib- 
rium solutions of isothermal self-gravitating clouds exist 
only above a temperature Tnim- For a sphere, Tmin is the 
temperature corresponding to the critical stable cloud, from 
Eq. inifor the fixed mass. For a cylinder, Tmin derives from 
Eq. |5l for the fixed mass line density. Mathematically we 
can express ratio of the temperature to Tmin in terms of 
/ and so the gravitational state of the clou d is given by 
r/Tmin > 1. For spheres T/Tmin = l/V/sph and for cylin- 
ders T/T„,i„ = l//cyl. 

The cloud parameters (diameter, central density, and 
column density) for both cylinders and spheres are shown 
in Fig. ini where as a minimum temperature we have chosen 
Tmin — 10 K. As the cloud cools at fixed mass line density 
and pcxt , its size shrinks while both central density and the 
column density through the centre rise. The variation of 
the radius and central density with temperature is stronger 
in the case of a cylinder compared to a sphere. For a given 
temperature ratio, the central column density through a 
cylinder is larger than that through a sphere if we neglect 
the complication of supercritical stable cases. 



* An interesting consequence is that if in a molecular cloud a 
large fraction of gas were bound in already near-critical stable 



cores, a pressure enhancement caused by a supernova shock or 
stellar wind would lead to a sudden burst of star formation in 
the whole cloud. 



8 



Jorg Fischera and Peter G. Martin: Physical properties of interstellar filaments 




10 1 10 1 

T IT =1 1 f T IT = f""^ 
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Fig. 9. Physical parameters of an isothermal self-gravitating cylinder (left panel) and sphere (right panel) with fixed 
Peyit/k = 2 X 10'' K cm~^ as a function of cloud temperature. Plotted are the cloud diameter (2rci, dotted line, scale on 
left axis), the central density (dashed line, right axis), and the column density through the cloud centre (solid line, right 
axis). The values given for the size and the density depend on the minimum temperature Tmin assumed to be 10 K. For 
spherical clouds the parameters are shown for all physical solutions including the supercritical ones characterized through 
an overpressure larger than Pc/Pext ~ 14.04. The grey dashed-dotted line in the left hand figure gives the temperature 
ratio for fcy\ = 0.75. 



As it cooled to Tmin, a cylinder would approach asymp- 
totically a singular state of an infinitely thin cylinder with 
infinite overpressure, assuming it was not in the meantime 
fragmented through instabilities. At a lower temperature 
still, there would not be an equilibrium solution at all: it 
would collapse to a spindle. 

Through cooling a sphere reaches the critical point be- 
yond which further cooling would cause a pressure drop at 
the cloud edge so that the cloud would collapse in a free fall 
time. The physical solutions presented in Fig.[3]for pressur- 
ized supercritical spheres correspond to cloud temperatures 
above the critical temperature; these supercritical spheres 
are not accessible simply by cooling at fixed mass and ex- 
ternal pressure. 

For fixed low mass or low mass line density Tmin might 
lie at a value that cannot be reached by any efficient cool- 
ing mechanism. For example Bok globu les show tempera- 
tures not much lower th an typical 10 K (jMvers et al.lll983t 
iBenson fc MversI Il989l) . See also Appendix [Ol Clouds 
with Tmin below this empirical physical limit would there- 
fore remain as gravitational stable equilibrium configura- 
tions. A filament might possibly enhance its mass line den- 
sity through accretion, raising T-adn and obviating this cool- 
ing barrier. Such spheres would also need to accrete more 
material to become physically unstable and/or they might 
become unstable through increased compression which 

rr 1/4 

raises T^in as . 



4. Observed physical parameters of filaments 

In the following we validate the model by comparing the 
theoretical prediction with recent observations. This is im- 
portant to confirm because then application of this model 
to interstellar filaments would allow an independent esti- 
mate of astrophysical quantities. For example, it could be 
used to estimate the interstellar pressure and the distance 
to cloud complexes and star formation regions, and it could 



be used in combination of measurements of the dust emis- 
sion to estimate the dust emissivity. 

Imaging submillimetre observations with Herschel pro- 
vide estimates of the fwhm of the filaments, their den- 
sity profiles, and of the central column densities. Some 
uncertainties relating to such measurements are discussed 
in Appendix [Xj Here, we discuss the results for individ- 
ual filaments in Pol a ris, IC 5146, and Aquila presented by 
lArzoumanian et "aD (|201lt) . To describe the filament pro- 
files, they adopted a generalization of Eq. [5]for the density 
profile, parameterized by index a: 

p(r)=p,(l + (r/V8ro)2)-"/2. (36) 

At large cloud radii (r :s> \/8ro) the profiles of the filaments 
in IC 5146 indicated a ^ 1.6±0.3, less steep than the value 
four they expected for clouds in a vacuu m, in agreement 
with p revious measurements discussed bv lFiege" fc Pudrit^ 
(j200Gt ). Our interpretation is simply that an equilibrium 
cylinder that is pressure confined is the relevant model. 
Relatively flat profiles are a natural outcome for cylin- 
ders with a mass ratio /cyi considerably smaller than unity 
(Section [S31). 

In Fig. [TOl the deconvolv ed fwhm and A'^h(O) results 
of lArzoumanian et al.l ()2011h are compared with the self- 
consistent prediction of our model of pressure-confined 
cylinders, including the asymptotic behaviours. The param- 
eters for the filaments of the three observed regions occupy 
different areas in the fwhm-NulO) plane. The filaments in 
Polaris show the lowest central column densities with most 
values below A'h(O) < 10^^ cm~^. The sizes are also ap- 
parently smaller compared to the two other samples. The 
filaments in IC 5146 have intermediate central column den- 
sities with a mean A'h ^ 6 x 10^^ cm^. The highest column 
densities were observed for filamentary structures in Aquila 
with values up to ~ 10^^ cm^^. Overall, the data indicate 
a trend that filaments with low column densities are on 
average smaller than those with high column densities. 
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Fig. 10. fwhm-Nii{0) relation of filaments. Observed val- 
ues in Polaris (dark grey symbols), IC 5146 (black sym- 
bols) , and Aquila (light grey symbols), from Fig. 7 of 
I Arzoumanian et "all ([2011). Note that we are plotting the 
column density of H nucleons rather than H2 . We also low- 
ered their values by 2.33/2.8 to correct for the value of 
p, that they adopted. The horizontal dashed-dotted lines 
give the imaging resolution achieved in the three fields 
(black: IC 5146, dark grey: Polaris, light grey: Aquila). 
Predictions from the model of pressure-confined isothermal 
self-gravitating cylinders are given for three different incli- 
nation angles (0°, 60°, and 80°; solid black and grey curves). 
The cloud temperature of the molecular gas is taken to be 
T = 10 K. The external pressure provided by the ambi- 
ent medium is assumed to be Pcxt/fc = 2 x 10* K cni^^ 
(see discussion in Appendix IA.3[) . The dashed lines are the 
asymptotes for low and high overpressure Pc/Pext or /cyi 
(this is a logarithmic version of Fig. |8]). For an overpres- 
sure corresponding to that of a critical stable sphere, and 
for inclination 0°, the expected value is shown as the filled 
black circle. Likewise, the vertical dotted line indicates the 
column density of a cylinder with mass ratio /cyi = 0.75; 
the upper non-linear scale indicates other values of /cyi for 
this inclination. The expected range for a pressure regime 
1.5 to 5 X 10* K cm^^ is shown by light grey shading. 



As we see in Fig. [10] the model of pressurized cylinders 
shows a good agreement with the observations, although 
obvious deviations do exist. The observations are broadly 
consistent with a cloud temperature T = 10 K and an am- 
bient pressure of Pext/fc = 2 x 10* K cm^'^. Such a pressure 
is in the range expected, as discussed in Appendix lA.3l The 
data show a dispersion around the theoretical curve which 
can be attributed to a variation of the external pressure 
alone. Most data points lie in a pressure regime PoKt/k from 
1.5 X 10* K cm~^ to 5 X 10* K cm~^. Note how the posi- 
tion of the maximum lies on a line parallel to the high /cyi 
asymptote. An additional horizontal deviation of data from 
the model is expected because of the distribution of incli- 
nation angles; as shown in Fig. [TUl higher inclination could 
in principle explain some of the high values of the column 
densities. Note that at constant pressure, the fwhm-Nii{0) 
curve translates vertically with T. 



iKandori et al.l (|2005[ ) studied the physical parameters 
of Bok Globules using an idealized model of a pressurized 
isothermal self-gravitating sphere to fit the column density 
profile obtained from extinction measurements. The derived 
external pressures ranged from 2.1 x 10* K cm~^ to 1.8 x 
10^ K cm~^ with a mean of 5.7 x 10* K cm"^ There is 
overlap with the pressures derived here but on average they 
are higher. 

4.1. Behaviour at small column density 

Comparison of the data to the asymptotic behaviour at low 
/cyl (Eq. IC.3I) constrains the ratio T/pcxt- In combination 
with this, the position of the maximum (Eq. I28p constrains 
T. At low column density, a higher pressure could be ac- 
commodated by raising T, but for higher column densities 
near the predicted maximum the fwhm might become too 
large. Likewise, if no efficient cooling mechanisms exist to 
cool the cloud temperature below 10 K, then we have a 
lower limit on the required external pressure close to the 
cited values. 

The parameters of most of the filaments in Polaris seem 
to follow the relation expected for rather low mass ratio, 
/cyl < 0.1. Those structures are therefore not strongly self- 
gravitating and so possibly transient density enhancements. 
If they exist, filaments with even lower /cyi would be small 
compared to the instrumental resolution and of low con- 
trast in the im ages. Stronger back ground fluctuations (cir- 
rus noise; e.g.. iMartin et aLll2010f ) in the other two cloud 
complexes would produce a more challenging limitation to 
identifying faint low column density structures, and so the 
fact that filaments with low mass ratio are less frequent or 
entirely absent in the two other samples might result from 
this selection effect. 



4.2. Behaviour at large column density 

We think that inclination effects are unlikely responsible for 
the extremely high column densities reported. High column 
densities without the expected small sizes might instead be 
related to cold filaments which are embedded and pressure- 
confined within larger filamentary structures, as mentioned 
in Sect. 13.6.21 Our simple idealized model is therefore not 
easily applicable. 

For example, in the IC 5146 sample, the fila- 
ment with the highest c olumn density (numbered 6 by 
lArzoumanian et "all 1201 it) also has the highest reported 
mass line density (152 Mq pc^^), much higher than the 
expected maximum value for a filament with K correspond- 
ing to 10 K. This is a very interesting but complex filament 
with large variations in central column density and fwhm 
along the sinuous ridge. The images and average line pro- 
file show one or more cold irregular inner filaments inside a 
broader structure. The narrow filaments are often paired in 
parallel segments. How high the mass line density is judged 
to be depends on how far out radially the column density 
is integrated and this depends on where the adopted model 
profile meets the "background." Measurement of the em- 
bedded filament (s) above the complex background of the 
broad filament would produce physical parameters closer 
to the simple model. A slightly higher K would help too, 
moving the predicted curve vertically in Fig. 1101 



10 



Jorg Fischera and Peter G. Martin: Physical properties of interstellar filaments 



Our sense is that this whole filament is not in free fall 
radial collapse, in which case the self-gravity of the larger 
embedding structure needs to be balanced too; an isother- 
mal model is probably not appropriate, but a higher K (or 
perhaps even extra magnetic pressure) would be needed. 
This is open to experimental investigation by molecular 
velocity and line- widt h measurement s such as reported 
by Pin eda et all ()2010D . iKramer et al] (|1999[ ) have in fact 
mapped a 0.4 pc square region on the ridge of this fila- 
ment in some lines of the rare isotopologues of CO. Even 
within this small region they find many clumps and two 
ridges with differing velocites and line widths. Overall this 
dense central region, while clumpy, is fairly cold in terms of 
K though perhaps warmer than 10 K. On the other hand, 
for ^''CO on a larger scale (2. 5 pc) the velocity di spersion 
cr w 1 km s"^ (position C3 in lDobashi et al.lll992D . While 
perhaps exaggerated by optical depth effects, this is sugges- 
tive of K being an order of magnitude larger in the broad 
embedding filament. The situation seems to be similar for 
the widest filament (numbered 12) where the profile sug- 
gests an interior cold (sometimes paired) filamentary struc- 
ture surrounded by a more turbulent gas. 

Interestingly, no filaments with high central extinction 
(A^h(O) >^ 2 X 10^^ cm^^) and the corresponding low size 
{fwhm < 0.08 pc) were found. There is no obvious selection 
effect against finding these high fcy\ filaments though, as 
suggested above, measurement of their parameters might 
be an issue depending on their environment. 

If strongly self-gravitating filaments, those above a cer- 
tain mass ratio, do not exist, then why? This might be 
inherent to the process responsible for the formation and 
evolution of filaments generally, or might be because of their 
disruption and/or short lifetimes. While interesting, the for- 
mer possibility is beyond the scope of this paper. One clue 
to the latter possibility is that according Fig. [TU] the em- 
pirical cut-off in /cyl is near 0.75, where the cylinders have 
about the same overpressure as a critical stable sphere and 
so if appropriately fragmented might form structures that 
would collapse. This is explored further in the following 
section. 



5. Structure formation along filaments 

In most cases observed filaments show substructure, con- 
tributing to the difficulty of measuring their parameters. 
Most interesting for understanding the potential role of fil- 
aments, as opposed to just high column and volume density, 
in the star formation process are the condensed structures 
extracted, "prestellar cores" having strong self-gravity. By 
definition they ar e cool, not having de tectable emission at 
70 jim or 24 /Ltm (jKonvves et al.|[20Tol ). and so while they 
have the potenti al for forming prot ostars or young stellar 
objects (YSOs) (| Andre et all 120101 ) they are not yet vig- 
or ously doing so. T hey are being externally heated (stage 
E, iRov et alJl20TTh with httle internal energy being gener- 
ated. These "prestellar cores" appear in the higher column 
density filam ents at typically Ay > 8 mag, as in Aquila 
and IC 5146 (| Andre et alH2010HArzoumanian et al.ll20ll . 
Cold prot ostars are rarer and also often associat ed with 
filaments (jAndre et al.ll2010l: iBontemps et al.ll2010[ ). 

Even filaments with much lower column density, those 
with low /cyl such as in Polaris, have substructure. These 
are dubbed "starless cores" as well, potentially confusing 
nomenclature: although they have sizes typical of what de- 



fines cores (radii less than a few tenths of a parsec) they 
comprise much less than a solar mass and are not signifi- 
cantly gravita tionally bound and so do not seem destined 
to form stars (jAndre et al.ll2010f ) . 



5.1. Equilibrium filaments 

If the basic structure is described by an equilibrium fila- 
ment, as we have argued, then the question arises as to 
the origin of the substructure. Equilibrium filaments are 
subject to gravitational instabilities that are axisymmet- 
ric, producing structure al ong the axial d irections. Two 
types may be distinguished (|Nagasawalll987h . For high /cyi 
(?'cyi ^ V^tq) the instability forms out of compression. A 
local density increase raises the self-gravity and lowers the 
gravitational energy, leading to a growth in the amplitude 
of the structure. For low /cyi, those with a fairly flat profile 
(''cyl <C -v/Sro), the instability relates to deformation of the 
surface; in such a "sausage instability" it is the volume en- 
ergy of the surface integral of the external pressure that is 
lowered. Transitional behaviour mixing these cases occurs 
near /cyi = 1/2 (rcyi = V^ro). 
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Fig. 11. Length and time scales relevant to fragmenta- 
tion as a function of /cyi, in units of the fwhm (left axis) 
and l/^/AnGpc the sound crossing time of tq (right), re- 
spectively. The shortest wavelength growing mode is Acr 
(dashed-multiple dotted lines) . The fastest growing mode is 
Am (dashed-single dotted lines) which, in the linear regime, 
develops on a time scale (long-dashed lines). For com- 
parison, Acr,sph is the length along the filament required to 
encompass a mass equal to that of a critical stable sphere 
for the same Pext and K (solid line) . The black broken lines 
correspond to non-magnetized cylinders. The related grey 
lines show the increasing effect of a magnetic field and corre- 
spond to — 1%, 10%, and 100% of equiparti tion. The val- 
ues are based on the work of lNagasawal ()1987[ ). see App. |El 
The vertical dotted line marks the mass ratio /cyi = 0.75. 



Linear analysis of the dispersion relation (iNaeasawal 

|1987[ ) shows that the critical wavelength above which per- 
turbations grow is, for /cyi 1, 



Ac 



= 3.96 VSro = 3.96rcyi-y/(l - /cyi) = 2.58 fwhm^^i 



(37) 
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and, for /cyi 0, 

Acr = 5.87 rcyi = 5.87 y/K7iV8rQ = 3.39 fwhm^y^, (38) 

with a smooth variation for values of /cyi in between. The 
wavelength corresponding to the maximum growth rate. 
Am, is about twice Acr (a factor 1.96 and 1.84 for the above 
limiting cases, respectively). Figure [TT] shows the /cyi de- 
pendence of Acr and Am, relative to the fwhm. At least 
initially, the fastest growing mode is quite elongated, with 
Xy^/fwhm^yi ~ 5. This would also be the separation of the 
cores, in units of fwhm. But in absolute terms, both Am 
and fwhm become arbitrarily small in the two limits. This 
is shown in Figure fT2l 
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Fig. 12. Like Fig. [TTl but with lengths and time scales 
expressed in physical units. The molecular cylinders are 
assumed to have temperatures of Td = 10 K and to be 
pressurized by a medium with Pcxt/k — 2 x 10^ K/cm"^. 
Also shown is the corresponding fwhm (see Fig. [7]). 

An important issue is the growth time scale, which is 
proportional to l/^/AirGpc = ro/^/K, the sound crossing 
time across rg. As shown in Fig. [TTl for the fastest growing 
mode the constant of proportionality falls from 4.08 to 2.95 
as /cyl increases from to 1 (Nagasawa 1987). It is also 
convenient to define a time scale TKp = ^/K/(Ai^GpZ^ , 
that evaluates to 0.4 My for the adopted values of K and 
Pext- Thus for /cyl -J> 1, 

Tm = 2.95(1 - /cyi)rKp. (39) 

This clearly arbitrarily short. For /cyi near the time scale 
is 

Tm = 4.08TKp. (40) 

While this is much faster than the crossing time rcyi/-\/K, 
because rcyi ^ r-Q, it is still finite, of order 1.6 My for the 
canonical values, even longer for a warmer or more turbu- 
lent gas. This would set up "sausage instabilities." The time 
scale Tm in My is shown in Figure [121 

The effect of a uniform magn etic field a l ong t he fil- 
ament has also been studied by iNagasawal (|1987[ ). The 
"compressive instability" (finite /cyi) is not much affected 
because the field is axial, but the "sausage instability" 
(low /cyl) is strongly suppressed because the field fimits 



the surface distortion. In that case, both Am (see Fig. (TTj) 
and the time scale Tm are increased significantly, even for 
a field much smaller than equipartition strength (where 
p2 = B'^/{A-KpcK) = 1). Thus, in the ISM this instability 
seems unlikely to be relevant to determining the substruc- 
ture seen in filaments with low mass fraction as in Polaris. 

The question of the final outcome of the non-linear 
growth of the "compressive instability" is of course critical. 
ICurrvl (I200flh found a sequence of equilibrium structures 
within filaments. Their shape is prolate along the cylinder 
axis, within a tidal lobe within which the equidensity con- 
tours are closed. These equilibria are gravitationally dom- 
inated and so are rounder as the pressure (density) con- 
trast between centre and lobe increases (limiting axial ratio 
about 1.7). The behaviour of the lobe radius and the mass 
contained within the lobe with increasing pressure contrast 
is reminiscent of the Bonnor-Ebert sequence, though it is 
not known whether these cores become unstable to collapse 
at similar overpressures. 

The non-linear gr owth of the in stability has been fol- 
lowed numerically by llnutsuka fc M ivama (1997) for two 
cases of equilibrium cylinders with /cyi = 0.9 and 0.2. For 
the former they find a complete disintegration of the fil- 
ament into well separated spherical clouds suggestive of 
gravitational collapse. For the latter a distinctive spheri- 
cal core forms as well, with the central density increased by 
a factor nine, but it does not collapse. 

5.1.1. Mass of the fragment 

The maximum mass available to be concentrated by a per- 
turbation of wavelength A is 

Mfrag = \M/l - 2V8/cyl(l - /cyl)-| , ■ (41) 

This is deliberately cast in the same form as Eq. [51 for the 
mass of a critical stable sphere, where the prefactor is sim- 
ply 4.191. 

Even if the growth of the perturbation is effective in 
concentrating material into a fairly spherical core, it seems 
intuitive to us that the outcome is unlikely to be gravita- 
tional collapse unless the mass of the core is comparable 
to the critical Bonnor-Ebert mass for the external pressure 
confining the original cylinder. To accumulate the critical 
mass, the required wavelength Asph,cr is, when /cyi is close 
to 1, 

Asph,cr//w/imcyi = 4.191/ (4^8^22/3 - l/cyi(l - /cyi)) (42) 
and, when /cyi is close to 0, 

Asph.cr//w/imcyi = 4.191/ (^2V3V8/cyiy^/cyi(l-/cyi)) .(43) 

Thus for both of these limiting cases the perturbations 
would initially have to be extremely elongated compared 
to the fwhm. The locus of Xsph.cr/ fwhm in Figure [TT] 
contrasts these extreme values with the finite value of 
X^/ fwhm^yi (in the range 5-6). Thus in both limits of 
/cyl, any condensations resulting from perturbations would 
be of low mass compared to the critical Bonnor-Ebert mass; 
this is not a promising scenario for forming prestellar or 
protostellar cores. 
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Indeed none are observed for filaments with, low /cyi. 
Also, because of suppression by the magnetic field, we are 
less interested in the case fcy\ close to 0. It is hard to imag- 
ine a perfectly uniform filament being created in a turbu- 
lent medium; whether the structure that does exist is just 
a natural by-product of formation of the filament needs to 
be investigated. 

For filaments with very high overpressure, or /cyi close 
to unity, the implication is that a series of small dense con- 
densations would develop, each of insufficient mass to col- 
lapse. On a longer time scale these might ni erge to form 
more massive cores (jlnutsuka fc Mivamalll997t) , perhaps ar- 
riv ing at t he hig h density equilibrium configuration found 
by ICurrvl (j200ClD . According to the findings presented in 
Fig. [To] this /cyi — ?► 1 limiting case has little relevance to 
the interstellar medium as the precursor narrow high col- 
umn density filaments have not been observed. How far the 
results are biased through the method applied to charac- 
terize highly opaque filaments need to investigate as well. 

However, for intermediate /cyi we see in Figure [TT] that 
Asph.cr can actually be less than Am (and even Acr) and com- 
parable to the fwhm. This is a more propitious situation, 
starting from a much less elongated perturbation, to ac- 
cumulate a mass comparable to the critical Bonnor-Ebert 
mass. 

For such intermediate /cyi, Tm « 3.5(1 — /cyi)TKp. This 
is sufficiently rapid (Fig. [l2]) that initially uniform fila- 
ments could develop structure on astronomically relevant 
timescales. As discussed above, it seems plausible that some 
of these could accumulate a critical Bonnor-Ebert mass 
and become the observed cold protostars. If this were the 
explanation, then the protostars ought to be separated 
by about Am, about five times the fwhm of the embed- 
ding/undisturbed filament (Fig. [TT]). This appears to be on 
the order of what is observed (see below) but needs to be 
carefully quantified. 

5.2. Pearls on a string 

According to iMen'shchikov et all ([201Clt ). the above- 
mentioned prestellar cores, observed along high column 
density filaments, are about the same size as the width of 
the filamen t. Inspection of the image s of source positions 
prese nted (jMen'shchikov et al.l 120101 : lArzoumanian et"al] 
1201 ih shows that they are also separated by about the 
fwhm of the filament. Therefore, remarkably, they are 
arranged "almost like pe arls on threads in a necklace" 
(jMen'shchikov et al.ll20l61 ). Because of the small separa- 
tion it appears unlikely that those structures are caused 
directly by the gravitational instabilities discussed above. 
We have seen that condensations resulting from that ori- 
gin would have a separation of at least 2.5 times the fwhm 
(for intermediate mass ratios /cyi >~ 0.5). We have ex- 
amined t he above-mentioned filam ent 6 in IC 5146 in both 
SCUBA toi Frances co et al.l l2008l ) and Herschel archival 
submillimetre images. There are many clumps that would 
be extracted as prestellar cores, but not all are equal, and 
our impression is that these are somewhat clustered with 
major concentrations typically separated by 2' to 5', cor- 
responding to 2.2 to 5.5 times the given fwhm. The YSOs 
associated wi t h this filament found in the Spitzer study by 
iHarvev et all ([2008D are also well separated. This is perhaps 
more in line with what would be expected from gravita- 



tional instabilities, but then a hierarchy and perhaps time 
sequence of fragmentation is suggested. 

At this point, we feel that the origin of the substructure 
is an open question. Other than gravitational instability, 
perhaps it is a natural consequence of the formation of such 
filaments; in the presence of strong self-gravity, it might be 
difficult to form a perfectly straight and uniform cylinder. 
In our experience it is also a challenge to extract and char- 
acterize "sources" found as closely spaced inhomogeneities 
along a filament. 

Because a separation of one fwhm is less than Acr.sph 
(Fig. [TT]) . these prestellar cores, while strongly self- 
gravitating, would probably be stable against collapse, con- 
sistent with their being starless. Whether these presently 
starless clouds will form any stars in the future and there- 
fore be relevant to the initial mass function of the stars is 
also an interesting open question. 

5.3. Filaments with high mass line density 

For a filamentary structure with more than the maximum 
mass line density (Eq. [T]), there is no equilibrium solution. 
If they could be created, somehow, because of their high 
density they would collapse on a rapid timescale. 

There are certainly filamentary configurations observed 
that have a large mass line density, considerably greater 
than the maximum line density if evaluated for K corre- 
sponding to a low value like 10 K which we find appro- 
priate to simple narrow filaments. lAndre et al.l ([2010) con- 
clude that the gravitational fragmentation of such "super- 
critical"[3 filaments is responsible for the formation of the 
observed self-gravitating prestellar cores and protostars. 
There are two problems with this scenario. 

First, our interpret a tion o f the numerical calculations 
of jlnutsuka fc Mivamal (|1997D is that there is fragmenta- 
tion only if the mass line density is finely tuned to being 
just above the maximum value and/or if the initial condi- 
tions already contain large perturbations, which begs the 
question. Perhaps a relevant configuration could be set up 
beginning with an equilibrium filament with /cyi close to 1 
and then making it "supercritical" by lowering the temper- 
ature. Otherwise, the predominant outcome is radial col- 
lapse of the cylinder to a spindle. Eventually the collaps- 
ing spindle would become optically thick and could develop 
substructure; but this would be on a small scale because of 
the large density and small size of the collapsed filament 
(Inutsuka fc Miyama 1997). 

Second, in the process of forming the fragments, a "su- 
percritical" filament itself would collapse, be consumed, and 
disappear on the same rapid time scale. However, the ob- 
servations indicate that prestellar cores, and even YSOs, 
occur along identifiable filaments of finite size, even if it 
appears that the mass line density is "supercritical." This 
coexistence is a challenge to be satisfied by any model of 
the structure formation, even for initially equilibrium fila- 
ments. 

While high column and volume density are arguably 
fundamental to any process of fragmentation, our conclu- 
sion is that protostar and cluster formation must be a more 
complex process than simply the gravitational fragmenta- 
tion of "supercritical" filaments. Based on the complex ra- 

^ Despite our comments in Sect. 12.21 their terminology is 
adopted here in this subsection. 
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dial profiles of tliese high mass line density filaments, we 
think that some progress might be made by considering 
cold embedded filaments within a broader structure sup- 
ported by a gas with higher effective in which case the 
overall configuration is not actually "supercritical" and in 
free fall. 



6. Summary and discussion 

We have analyzed the physical properties of interstellar fil- 
aments on the basis of an idealized model of isothermal 
self-gravitating infinitely long cylinders which are pressur- 
ized by the ambient medium. The pressure-confined cylin- 
ders have a mass-line density smaller than cylinders ex- 
tending into a vacuum. The mass fraction /cyi, the ratio of 
the mass-line density to the maximum possible for equilib- 
rium configurations, describes the gravitational state of the 
cylinder. For given temperature and external pressure, it is 
used in deriving analytical expressions for the central den- 
sity, the radius, the column density profile, the central and 
average column densities and the fwhm. The dependence of 
the physical properties on external pressure and tempera- 
ture is clear in these expressions. The results are compared 
to the case of pressure-confined isothermal self-gravitating 
(Bonnor-Ebert) spheres, characterized by the mass fraction 
/sph, the ratio of the mass to the mass of the critical stable 
sphere. 

We compared the model prediction for the relation be- 
tween the size and the central column density with recent 
observations of filaments. Given the complexity seen it is 
gratifying to find good agreement, even surprising consid- 
ering the idealization of the model. In practice the model 
seems best applicable to those filaments that are cold and 
show a simple smooth elongated structure. On the basis 
of the model the filaments appear to experience an am- 
bient pressure Pcxt/fc in the range 1.5 x lO'* K cm^'^ to 
5 X 10^ K cm~'^ and to have a gas temperature T w 10 K. 
For these parameters, the observations indicate an apparent 
physical upper limit to the mass ratio /cyi ~ 0.75; higher 
column density narrower filaments are not seen. The cor- 
responding maximum central extinction for the observed 
filaments fitting the model is Ay « 7 mag. This is interest- 
ingly close to the empirical threshold for star formation, but 
we think this is coincidental, es pecially given the variable 
threshold from region to region (jEnoch et al.ll200"7l ). 

We have summarized and discussed previous stability 
studies that have shown that all unmagnetized infinitely 
long cylinders /cyi < 1 are subject to gravitational instabil- 
ities. While a magnetic field can have a stabilizing effect 
for /cyi < 0.5, in particular suppressing "sausage insta- 
bilities" for /cyi « 0, there is a negligible effect on com- 
pressive instabilities for /cyi > 0.5. Interstellar magnetized 
filaments would seem to become increasingly unstable to- 
wards large /cyi because of a strongly increasing growth rate 
of the instabilities. The disintegration of filaments close to 
the maximum mass line density (/cyi > 0.9) would produce 
fragments whose masses are well below the critical value 
to form stars. We found no indication in the theory for a 
strong threshold for the formation of fragments that would 
collapse. However, intermediate values of /cyi seem most 
favorable. The fastest growing disturbance for /cyi > 0.5 
occurs at a length approximately five times the fwhm, so 
that protostars or clusters would be well separated. 



It is difficult to understand in this model why in high 
column density filaments there would be substructure akin 
to the extracted prestellar cores separated along the fil- 
ament by only one fwhm, which is well below the mini- 
mum length at which a disturbance would grow (Acrit >~ 
2.5 X fwhm^y]). This might point to external influences ac- 
companying the process that formed the filament in the first 
place, resulting in the non-uniform mass line density along 
these filaments, not the idealized model considered here. On 
the other hand, if one concentrates on the major mass con- 
centrations along the filament, their separation might be 
more compatible with an origin in gravitational instability 
along the cylinder. 

High mass line density filaments are observed that 
would appear to have more than the maximum mass line 
density that could exist in equilibrium, at least if evalu- 
ated for the low adopted K. The coexistence of embedded 
narrow high column density filaments, even in pairs, and 
considerable substructure (prestellar cores) and protostars 
suggests to us a scenario more complex than simply the 
free-fall collapse and fragmentation of cylinders with more 
than the maximum line density. Although the simple model 
cannot explain such complex structure, it points to a sit- 
uation in which the embedding structure is supported by 
gas with a K larger than the low value that might reason- 
ably be adopted for the narrow embedded structure. This 
is obviously open to observational scrutiny via mapping of 
velocit ies and li n e wid ths of appropriate molecules. 

As iMcCreal (|1957D has emphasized, "genuine" gravi- 
tational collapse to form stars and clusters is best real- 
ized in configurations with roughly the same dimensions 
in all directions so that gravitational effects are three- 
dimensional. He argued that the breakup of less favorable 
one-dimensional configurations like filaments must be due 
to irregularities in the density distribution and the exter- 
nal pressure to which they are subjected, or to differential 
motions. Empirically, there is suggestive evidence that such 
breaking of the symmetry of an idealized long cylinder is 
important, as occurs for example where filaments appear 
to cross. These are matters that could be explored now in 
the results of numerical simulations. 
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Appendix A: Questions of scale 

In comparing the models with the data we need to be aware of the 
context, what can be measured and the related uncertainties. In this 
appendix we discuss column density, size, external pressure, and ef- 
fective temperature. 

A.l. Column density 

Column densities of filaments from Herschel data are fundamentally 
in units of the dust optical depth in the submillimetre. This is con- 
verted to A'h using a submillimetre dust opacity. This opacity ap- 
pears to show a depe nde nce on environment and might be uncertain 
by a factor two (Martin et al. 2011). For the overview comparison in 
Figure [TOl where Nn{0) is on a logarithmic scale, this is might not be 
a major concern. 

In the diffuse and largely atomic interstellar medium, where both 
A'^H S'lid optical interstellar extinction curves can be obse rved directly, 
we ha ve the ratio Nii/E{B-V) = 5.8 X lO^^cm-^/niag l IBohlin et "akl 
1197^ and an absolut e to relative extinction Ry = Ay/E{B — V) = 
3.1 llFitzpatricklll999l ). A column density of lO'^'^ cm~^ therefore cor- 
responds to Ay = 5.3 mag. We use this same conversion when Ay is 
being used a "shorthand" for column density in dense molecular en- 
vironments, recognizing that the actual visual extinction, if it could 
be measured, would probably not be this value. 

Another measure more applicable at higher column densities is 
near infrared colour excess. This is also often converted to Ay, which 



again should be regarded as a "shorthand" because the shape of the 
interstellar extinction curve through the optical probably varies. It 
might be more reliable to convert the near infrared colour excess di- 
rectly to A^H, but that calibration is not actually know n for the high 
colum n densities encountered in star-forming regions l IMartin et al.l 
120111). 

Another uncertainty might be related to radiative transfer which 
produces a radial variation of the dust temperature within externally- 
heated filaments. While the dust opacity seems to be larger for dust in 
denser environments, the averaging of the dust emission over a distri- 
bution of temperatures leads to a reduced effective emission coefficient 
and a spectral dependence of the emissivity that is less steep than is 
intrinsic for dust in a diffus e environment as theoretically demon- 
strated for condensed cores llFischeral l2011f) . The radiative transfer 
effects increase with central column density and could therefore play 
a larger role at high mass ratios. The column density profile indicated 
by the surface brightness would be flatter than the intrinsic profile. 



A.2. Size 

Sizes of fllaments are measured from the proflle of the column den- 
sity in the direction orthogonal to the axis of the fllament. The fwhm 
could be measured directly but is more usually found as a parame- 
ter of a model flt to the proflle. The model flt also needs to account 
for the background on which the filament appears superimposed. In 
practice both the background a nd the profile ca n vary along the fila- 
ment, adding to the challenges l lArzoumanian et al.ll2011l resorted to 
using a Gaussian model to find the fwhm) and the reported standard 
deviation in the parameters. Even in nearby molecular clouds, the 
filaments, while resolved, are often not a lot larger than the imaging 
resolution. Therefore, the reported sizes have been "deconvolved." In 
the future, this might be avoided by forward modeling, adopting an- 
alytical models like we have presented, convolving the predictions to 
the resolution of the data, and then iteratively fitting the parameters. 

The background is still something to be flt too. Depending on 
the circumstances, this background could be from foreground or back- 
ground material unrelated to the fllament, and/or a contribution from 
the medium providing the pressure confinement. A well-dcflned edge 
of a fllament might suggest, for example, that the fllament can be 
regarded as cold material which is surrounded by a warmer or more 
turbulent but less dense environment. This relates to the following 
discussion of external pressure. 



A. 3. External pressure 

IBoulares & Coj3 (119901 ) presented a model of the vertical support of 
the Galactic disk. The mid-plane pressure near the sun is at least 
Poxt/fc = 2.5 X 10* K cm"''. About a third of this is the effective verti- 
cal support by cosmic rays; however, on the size scale of the structures 
of in terest her e, cosmic ra y pressure ca n be ignor ed I Curry fc McKe3 
HOOO); Fische ra & Dopita (2008); Fischeri (120111 '). The remainder is 
a combination of kinetic pressure (thermal and turbulent motions) 
and magnetic pressure. For a magnetic flcld of less than equipartition 
strength, the kinetic pressure would be at least 0.8 X 10-* K cm 
The thermal pressure is only about a third of this llKalberla & Kerpl 
120091 ). the rest being turbulent pressure. 

The filaments being discussed are within clouds with sufH- 
cie nt column density to have become molecular {{Ay) > 1; e.g.. 
Glo ver & Mac Low! 120111 ). If it is assumed that these clouds can 
themselves be described as pressure-confined isothermal spheres, with 
enough self-gravity to be near the critical mass, then following Eg. 1211 
and Eq. lB.lOl the external pressure required to produce a mass surface 
density (E) averaged over the projected radius is given by0 

Pcxt = (7r/8)G(S>2 or 
p^^t/k = 0.37 X lO"* {Ayf K cm-3. (A.l) 



Using a dimensio nal analysis of a cloud in virial equilibrium, 
iMcKee fc TanI (|2003l ') derive a similar relationship between the 
average internal pressure and the column density. There is a 
similar coefficient of proportionality, not surprisingly because 
both configurations are strongly self-gravitating and because the 
average pressure in a c ritical sphere is only 2.5 times larger than 
the external pressure (|Spitzedll96l ). 
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A subcritical spliere would require a greater external pressure to pro- 
duce the same average column density (Fig. The average col - 
umn density in the GIvIC sample discussed by iHever et al.l ll2009l i 
corresponds to (Ay) 2.8, which would require a pressure of 
2.9 X 10^ K cm~'^. This seems consistent with the fact that this cloud 
sample is in the inner Galaxy where the midplanc pressure is probably 
higher. 

The central extinction through a critical stable isothermal sphere 
is Ay = 8.0\/(pcxt/fc)/(2 X 10" K cm"^) mag (Fischera fc Dopital 
l2008t Fig. This predicted amount of extinction seems 

to be in agre ement with that observed for molecular clouds 
lIDopita fc Suth erland 2003). 

Two conclusions arc that the midplane pressure is indeed quite 
high, and that the pressure within the molecular cloud due to its self 
gravity — the ambient external pressure experienced by the filaments 
— will be higher still, depending on where the filaments lie within the 
pressure profile of the molecular cloud. 

Another contribution could be due to ram pressure, for example 
related to the sweeping up of the North Celestial Loop of which the 
Polaris field is a part. 



A. 4. Effective temperature 

For the model we need the effective temperature embodied in K, re- 
lated to both the thermal and the turbulent (non-thermal) motion of 
the gas. The most direct way to measure this would be through the 
line width of an appropriate mol ecular tracer. For exa mple, the ammo- 
nia mapping of B5 in Perseus by I Pineda et al.l ||201CI| ) revealed a sharp 
decrease in line width in the central 0.1 pc. This is thought to be the 
result of the decay of supersonic turbulence present on larger scales; 
the observed li ne widths ar e subs onic for a thermal/kinetic temper- 
ature of 10 K. [Foster et al.l i200gf) present results from an ammonia 
survey of many other cores in Perseus and find them mostly quiescent 
("■nonthermal ~ 0.12 km , but See their histograms for details of 
the full survey). The modeling of two lines gives kinetic tem peratures 
about 12 K; derived excitation temperatures are lower CFo ster et al.l 
f2009). The dust temperature from the submillimetre spectral en- 
ergy distribution observable by Herschel is typically low (down to 
10 K) in the high column density regions in low mass star form ing 
regions (see the temperat ure maps in, e.g.. IBontemps et al.|[20Tci and 
lArzoumanian et al.ll201l] i. Coupling of the kinetic temperature in the 
gas to the d ust t emperature comes into play for > 10* cm~^ 

ll Gold smith 2001). According to Fig. \T\ this becomes relevant for 
strongly self-gravitating cylinders (/cyi > 0.7). Taken together, while 
perhaps a platonic ideal, a small dense isothermal core with K corre- 
sponding to 10 K is reasonable as a reference value. However, K might 
be higher than this in some environments or stages of evolution. And 
of course the structure might not be exactly isothermal. 



for the cloud radius we obtain for the mean mass surface density of 
an isothermal self-gravitating sphere 



(2> (^ci) ^ 



where 



C(Zcl) 



A-fsph(^cl) _ V47rGpext 
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C(2cl), 
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We want to analyze the function (^(z) for the values z; with 
i = 1,2,... at which a spherical cloud of a fixed mass produces 
pressure maxima at the cloud outskirts for changing cloud size z^i 
llFischera &: Dopital l200"^ i. The global pressure maximum at zi = 
2max with Zmux ^ 6.451 characterizes the critical stable cloud which 
has an overpressure e'^^^™"'') ^ 14.04. The last two terms in Eq. IB. II 
for a critical stable cloud produce a value of approximately 4.191. 
We obtain the condition for the pressure maxima from Eq. lB.ll which 
we transform to an expression for the pressure at the cloud outskirts 
p(2cl)- At the pressure maxima we find that the derivative of the po- 
tential needs to fulfill the condition 



dz 



Jo 

On the other hand, from the Lane-Emdcn equation it follows that 



(B.6) 



dl 



dz'z'^e-'^f^'). 



(B.7) 



Using this expression to replace the derivative in Eq. IB. 61 we obtain 
for the condition for the pressure maxima 



-Lo{zi)/2 



(B.I 







For clouds producing maxima at the cloud outskirts the two last terms 
in Eq. IB.ll of the cloud mass are therefore equal to •/2 z? e~"(^i). For 
those clouds we get the alternative mass equation 



(B.9) 



For the critical stable cloud we have the condition 



~ 4.191 with Zmax ^ 6.451. 
If we insert the condition IB. 81 in Eq. IB. 51 we see that the pres- 
sure maxima are at z- values where C{zi) = \/2 and we obtain as the 
asymptote for large overpressure: 



Appendix B: Mean column density of spheres 

In paper II the approximation for the mean column densities of cylin- 
ders and spheres as a function of the overpressure were presented. 
Here, we want to comment on the values for spheres and give an 
exact expression for the asymptotic behaviour at large overpressures. 

The mass of a spherical isothermal self-gravitating cloud which is 
in pressure equilibrium with the surrounding medium (p(Z(.i) = Pcxt) 
is given by (paper I) 



sph 



.e-'"^/^ / dz'z'^e-'^t-'). 

^ 47rG3pext Jo 



(B.l) 



Here, w = <f>/ K = Inpc/Poxt is the unit free potential and z = rA the 
unit free radius of the Lane-Emden equation 



(B.IO) 



The asymptotic value of the mean column density given in Eq. 1211 is 
obtained using the relation 

(ATh) = (S> (AmH)-^ (B.ll) 

Appendix C: The fwhm of isothermal 
self-gravitating clouds 

To derive the fwhm we replaced the density through p(r) = pce~" 
and the radius by r = z/A. In case of cylinders the potential is given 
by aj(z) = 21n{l-|-z2/8} (Eq.[2ll. The mass surface density at impact 
parameter z < z^i is then given by 



1 d / 2 do; 
z^ dz V dz 



(B.2) 



S(z) = 2,/^e"(-ci)/2 / dz— L=f 
^ ' \l AtxG / Vzi - z^ 



(C.l) 



where A^ = 1/)"q = A-TzGpc/K. Using the expression for the mass 
(Eq. lBTTI l and 



''cl = 2clM = Zc 



K 



(B.3) 



The impact parameter Zf^i^^ for the fwhm for an isothermal pressur- 
ized cloud is determined by 



2 / dz- 



-uj(z) 



9 2 

z — Zf . 



dze-'^(^). 



(C.2) 
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For z^i <g; 1 the density profile is essentially fiat and we obtain 
the simple asymptotic behaviour zj^f^^ = y/s/iz^i which gives 



If we identify x as the value where the function y reaches its maximum 
we have 



fwhm(z,i <1) = ^ S(0). 

2 Pcxt 



(C.3) 



to 



b + 1 
a- 1' 



At sizes 2ei 3> 1 the right hand side approaches asymptotically a 
constant value. For spheres we get 



/oo = / dze-'^(^) ^ 3.028. 
/o 

For cylinders we obtain 
\z ' 



(1 + ^2/8)2 V2' 



(C.4) 



(C.5) 



We find that the combined function can be written as 

it/tor-' 



Vit) = S/maxC 



with 



(a - l)''-l(fe + l)''+l 



(a + b) 



a + b 
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(D.8) 
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This constant determines the impact parameter for high overpres- 
sure or large z^i. For spheres we find 



^fwhm ^ 2.997 

and for cylinders 



= \/8(22/3 -1). 
The relation r = z/A provides 
K 



fwhm{Zc\ > 1) = 2zf^hm 



V47rGpcxt 
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In the special case a = 2, 6 = as for the relation between fwhm and 
S(0) for cylinders or in good approximation for a sphere we have 

(D.IO) 

(D.ll) 

(D.12) 

(D.13) 




As noted in Sect. 13.51 in the limit of high overpressure {z^i 2> 1) the 
fwhm for constant T and poxt decreases inversely proportional to the 
square root of the overpressue pc/pcxt. Replacing the central pressure 
by the mass surface density through the centre gives finally for z^i ^ 1 
the asymptotic behaviour 



fwhm{zci > 1) = itlhULj^^ 



K 1 

S(0)' 



Appendix E: Analytic approximations of the length 
and time scales of perturbations 

The derived dependence of the perturbation lengths Acr and Am and 
the initial growth time Tm on t he mass rat i o and the magnetic field 
strength are based on the work of lNagasawal lll987f) who presented cal- 
culations for /cyi = 0.01, 0.5, and 1.0. We assumed that the mass ratio 
dependence of the time scale given in units 1 / ^/i-KGpc and the ratio of 
the disturbance length to the fwhm of the cylinder for non-magnetized 
cylinders can be described through a fourth-order polynomial function 



(C.9) 



y{x) : 



so that fwhm oc KT, ^(0). For high overpressure the relation fwhm- 
S(0) does not depend on the external pressure. 



Appendix D: On combining power laws 



i=0 



(E.l) 



Here, we discuss the main properties of the function used to fit the 
relation of the fwhm and the central column density A^hCO). Call these 
y and x, respectively. There are asymptotic limits relating these for fwhm^yi/ (VSro) 
small and large mass ratios (e.g., the power laws in Ea. l31l and Eg. 1291 
for cylinders). We assume that two power laws 



in a; = /^yi with vanishing derivatives for low and high mass ratios 
= 2/'(l) = 0)- The curves reproduce the theoretical values at 
/cyl = Oi 0-5i a^nd 1.0. The disturbance lengths can be transferred to 
absolute lengths using a fifth-order polynomial approximation for the 
fwhm in units of the characteristic length ^/Sro given by 



(E.2) 



yi 
y2 



X <^ X, 



(D.l) 
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with a > 1 and fe > — 1, can be combined in an interpolating formula 



y{t) = iis" 



(i + (tAo)^)°^'' 



(D.3) 



The approximation is 1% accurate and provides the correct asymp- 
totes at low and high /cyi- The coefficients ai for the polynomial 
representations are given in Table [ETT] 

To estimate the dependence of the magnetic field on the parame- 
ters Acr, Am, and Tm we used at low mass ratios the dispersion relation 
for an incompressible self-gravitating cylinder which has been shown 
to be a valid re presentation for f^y\ — > 0. The relation is given by 
l lNaeasawalll987f) 



where we have introduced t = x/x and where tg = {£,2/^,1)^^^°''^^^ /x. 
As can be verified, the equation provides the correct limits for t <^ 
to and t ^ to- The parameter 7 describes the smoothness of the 
transition and is adjusted to optimize the fit. For 7 — > 00 the function 
becomes a simple broken power law. The maximum lies at 



47rGpc /o(5) 



1 



8/cyl/0(O^l«) 



(E.3) 



where Iv{£) and Ki,{£) are modified Bessel functions, p^ = 
B^/iAwpcK), 5 = (27r/A)i?, and uj = l/r. The behaviour at larger 
/cyl we described through a fourth-order polynomial function where 



we have chosen x 



/^Yf for "Tm and X 



f. 



1/2 

cyl 



for the disturbance 



to 



1 



b+1' 



with 



a + b T 



(a - l)"-i(f)-f- 1)'' 
(a + 6)"+'' 



(D.4) 



(D.5) 



lengths Acr and Am. To the analytical function at low mass ratios we 
produced a smooth transition and assumed a vanishing de rivative at 
.fcvi = 1 ■ The curves were chosen to reproduce the values of lNaeasawal 
ll 19871 1 at /cyl = 0.5 and /cyi = 1. For strong magnetic fields the left 
boundary for the polynomial approximation is set to a value were we 
obtain at high mass ratios essentially a flat curve. The correspond- 
ing curves only show qualitatively at which /cyi the parameters are 
affected by the magnetic field. 
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Table E.l. Constants of polynomial approximations 



VSro 



ao ai 

4.08 0.00 

3.39 0.00 

6.25 0.00 

0.00 1.732 



a2 as 

-2.990 1.460 

-2.414 1.588 

-6.890 9.180 

0.000 -0.041 



a4 as 
0.400 0.00 
0.016 0.00 

-3.440 0.00 

0.818 -0.976 
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